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ABSTRACT 


The  serpulid  polychaete,  Hydroides  elegans  is  a  problem  fouler  in  warm  water  ports  throughout 
the  world  and  is  an  excellent  candidate  for  investigating  molecular  reception  during  larval 
settlement  and  metamorphosis.  H.  elegans  is  induced  to  settle  and  metamorphose  when  given 
the  appropriate  bacterial  cue.  This  study  utilized  transcriptome  sequencing  of  H.  elegans  and 
genome  sequencing  of  two  inductive  bacterial  strains  (i.e.,  Pseudoalteromonas  luteoviolacea  and 
Cellulophaga  lytica )  to  provide  a  means  to  examine  mechanisms  and/or  genes  that  may  be 
significant  in  understanding  the  interaction  between  bacteria  and  larvae  that  result  in  biofouling. 
Transcriptome  analysis  and  in-situ  hybridization  techniques  determined  the  genetic  expression 
and  location  of  genes  in  II.  elegans  as  well  as  important  biological  pathways  involved  in 
metamorphic  competence.  Likewise,  the  genome  sequencing  of  P.  luteoviolacea  and  C.  lytica 
determined  significant  coding  sequences  for  structures  that  play  a  role  in  the  settlement  and 
metamorphosis  of  II.  elegans. 
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I.  Scientific  and  Technical  Objectives 


The  tubeworm  Hydroides  elegans,  is  a  significant  problem  fouler  in  warm  water  ports  around  the 
world  and  is  an  excellent  candidate  for  investigating  molecular  reception  during  larval  settlement 
and  metamorphosis  (Hadfield  et  al.  1994).  Modern  molecular  techniques  enable  us  to  examine 
and  identify  larval  receptor  systems  providing  a  significant  contribution  toward  understanding 
the  interaction  between  bacterial  ligands  and  larval  receptors  that  result  in  biofouling. 

Expanding  upon  the  research  funded  by  ONR  Grant  No.  N000 14-08- 1-04 13,  we  will  determine 
the  molecular  basis  of  bacterial  induction  of  larval  metamorphosis  and  identify  biological 
pathways  and  genes  expressed  at  the  stage  of  metamorphic  competence  of  the  larvae  of  //. 
elegans.  Additionally,  H.  elegans  is  induced  to  settle  and  metamorphose  as  a  response  to  a 
bacterial  cue  from  a  bacterial  biofilm.  Recently,  this  cue  has  been  identified  to  be  a  phage-tail 
like  bacteriocin  produced  by  the  bacterium  Pseudoalteromonas  luteoviolacea  (Shikuma  et  al., 
2014).  However,  in  addition  to  P.  luteoviolacea,  Cellulphaga  lytica  has  also  been  identified  as  a 
bacterium  capable  of  inducing  the  settlement  and  metamorphosis  of  H.  elegans.  Therefore, 
sequencing  the  genome  of  these  inductive  bacterial  species  allows  discovery  of  bacterial 
products  and/or  metabolites  that  may  be  involved  in  the  induction  of  larval  settlement  and 
metamorphosis  -  the  fundamental  basis  of  marine  biofouling. 


II.  Approach 

A.  Transcriptome  Sequencing  of  Hydroides  elegans 

The  initial  approach  to  screening  H.  elegans  for  candidate  receptor  genes  involved  a  large 
arrayed  cDNA  printed  library  developed  in  the  laboratory  of  Dr.  Andrew  Cameron  at  the 
California  Institute  of  Technology.  Approximately  200,000  cDNA  gene  fragments  that  were 
cloned  in  bacteria  were  contained  on  a  filter  to  be  screened  for  possible  receptor  genes.  In  this 
process,  degenerate  primers  were  created  from  sequences  of  candidate  genes  that  were  not 
necessarily  specific  to  H.  elegans  but  instead  found  in  the  GenBank  database  of  other  marine 
invertebrates.  Recently,  however,  more  advanced  and  accurate  methods  for  “mining”  for  genes 
of  interest  have  developed  through  Next  Generation  Sequencing  (NGS)  technology.  One  such 
method  used  in  NGS  is  the  Illumina  HiSeq  method.  Illumina  HiSeq  allows  for  the  generation  of 
multiple  genetic  sequences  (i.e.,  reads)  from  a  biological  sample.  As  a  result,  the  subsequent 
approach  taken  by  this  laboratory  to  determine  the  potential  receptor  sequences  and  expression 
of  specific  proteins  involved  in  the  settlement  and  metamorphosis  of  H.  elegans,  was  to  sequence 
the  transcriptome  of  H.  elegans.  The  transcriptome  of  the  larval  (precompetcnt  and  competent) 
and  adult  stages  of  H.  elegans  was  determined  since  the  transcriptome  of  any  organism  reflects 
the  genes  that  are  being  actively  expressed  at  any  given  time. 

Specifically,  total  RNA  samples  were  prepared  using  the  Qiagen  RNeasy  Mini  Kit  and  the  RNA 
quality  was  determined  with  the  Agilent  Bioanalyzer.  The  total  RNA  sample  was  then  submitted 
to  the  Huntsman  Cancer  Institute  at  the  University  of  Utah  and  the  subsequent  isolation  of 
mRNA  was  used  for  Illumina  HiSeq  101  paired  end  sequencing.  Sequencing  for  all  larval  and 
adult  stages  was  run  in  triplicate  (biological  replicate)  and  in  individual  lanes  (technical 
replicates)  of  the  flow  cell.  The  transcriptome  sequences  generated  with  the  Illumina  HiSeq 


method  were  aligned,  annotated  and  assembled  with  specific  software  programs.  Additionally, 
those  genes  that  were  differentially  expressed  in  precompetent,  competent  and  adult  stages  were 
determined  through  bioinformatic  analysis.  Primers  were  designed  to  the  transcriptome 
sequences  to  amplify  potential  receptor  genes  and  then  used  in  the  following  capacities:  1)  to 
produce  a  700-900  basepair  fragment  and  2)  to  use  the  PCR  product  in  in-situ  hybridization 
studies.  The  riboprobes  that  are  generated  for  the  in-situ  hybridization  studies  were  then  used  to 
assess  the  anatomical  location(s)  of  the  transcripts  in  the  competent  larvae  of  H.  elegans  using 
the  Roche  Dig/a-Dig  labeling  system  and  visualized  with  a  fluorescent  microscope. 

Additionally,  the  genes  that  were  differentially  expressed  and  annotated  in  competent  larvae 
were  then  imported  into  the  program  DAVID  (Database  for  Annotation,  Visualization  and 
Integrated  Discovery)  to  determine  the  biological  pathways  and  genes  involved  at  the  stage  of 
metamorphic  competence. 

B.  Genomic  Sequencing  of  Inductive  Bacterial  Strains 

A  second  approach  involved  the  genomic  sequencing  of  two  bacterial  genera  that  are  capable  of 
inducing  the  settlement  and  metamorphosis  of  H.  elegans.  The  first,  more  widely  studied  marine 
bacterium  is  Pseudoalteromonas  luteoviolacea  (HI1)  and  the  second  is  Cellulophaga  lytica 
(HI1). 

Both  P.  luteoviolacea  (HI1)  and  C.  lytica  (HI  1 )  were  originally  isolated  from  a  seawater  table  at 
the  Kewalo  Marine  Laboratory  (Honolulu,  HI).  The  DNA  was  extracted  from  overnight  cultures 
using  the  MoBio  UltraClcan  Microbial  DNA  Isolation  Kit  and  DNA  quality  and  quantity  were 
determined  using  the  Nanodrop  2000  platform  (Thermo  Scientific).  The  isolated  DNA  was  then 
submitted  to  the  National  Center  for  Genome  Resources  (NCGR)  for  PacBio  Single  Molecule 
Real  Time  (SMRT)  Sequencing,  a  third  generation  sequencing  method. 


III.  Concise  Accomplishments 

A.  Transcriptomes  of  Hydroides  elegans 

Approximately  100  million  reads  were  generated  by  the  Illumina  HiSeq  101  paired  end 
sequencing  method  for  each  developmental  stage  of  H.  elegans.  The  reads  were  then  aligned 
using  the  software  program  Burrows  Wheeler  Aligner  (BWA).  Once  the  reads  were  aligned  with 
BWA,  the  de  novo  reconstruction  of  the  transcriptomes  from  the  RNA-Seq  data  was 
accomplished  by  the  Trinity  software  program.  Each  sequence  was  then  submitted  to  a  Blast 
search  against  various  repositories  of  protein  databases  found  at  NCBI,  Swiss-Prot,  TrEMBL  and 
NR  (Non-Redundant  database).  Also,  those  sequences  that  were  unique  to  each  stage  were 
determined  using  EBSeq,  a  software  program  for  differential  expression  analysis  and  categorized 
into  different  patterns  (i.e.,  Pattern  1-5)  (Table  1). 

The  initial  use  of  the  transcriptome  sequencing  data  was  to  design  riboprobes  that  targeted  genes 
related  to  the  innate  immune  system  since  Ii  elegans  responds  to  bacterial  cues  and  the  innate 
immune  reaction  is  the  first  line  of  defense  for  animals  against  various  pathogens  (i.e.,  bacteria, 
viruses,  etc.).  Thus,  the  riboprobes  were  designed  targeting  the  innate  immune  receptors; 


Galaetose  Binding  Leetin  (GBL),  Peptidoglyean  Recognition  Protein  (PGRP)  and  NF  Kappa  B. 
GBL  is  a  C-type  leetin  reeeptor  that  reeognizes  galaetose  in  bacterial  LPS.  PGRP  is  a  highly 
conserved  pattern  recognition  receptor  that  binds  to  components  of  the  bacterial  cell  wall  and  can 
activate  toll  reeeptors  as  part  of  the  innate  immune  response.  NF-kappa-B  is  a  protein  complex 
that  controls  the  transcription  of  DNA  (transcription  factor)  in  response  to  stimuli  such  as  stress 
and/or  viral  or  bacterial  antigens.  The  riboprobes  were  used  in  in-situ  hybridization  experiments 
to  locate  and  visualize  the  anatomieal  location  of  the  transcripts  on  the  competent  larvae  of  II. 
elegans. 

Additionally,  as  mentioned  previously,  the  genes  that  were  differentially  expressed  in  competent 
larvae  were  determined  (i.e.,  Pattern  3)  (Table  1)  with  the  software  EBSeq  and  analyzed  using 
DAVID.  Of  the  6,422  eoding  sequences  expressed  in  competent  larvae,  4,679  contigs  were 
associated  with  peptides  of  known  proteins.  Of  the  4,679  peptides,  203 1  peptides  were  unique  to 
competent  larvae.  DAVID  analysis  was  able  to  further  classify  the  2,031  competent  larvae  gene 
expression  patterns  into  KEGG  Pathways  and  Gene  Ontologies  (GO).  The  KEGG  Pathway  was 
designed  to  link  genes  in  the  genome  to  gene  produets  in  a  pathway  to  see  which  pathways  and 
associated  functions  are  likely  to  be  encoded  in  the  genome  of  the  organism.  KEGG  pathways 
are  classified  into  sections  whieh  include  1)  Metabolism  2)  Genetic  Information  Processing  3) 
Environmental  Information  Processing  4)  Cellular  Processes  5)  Organismal  systems  6)  Human 
Diseases  and  7)  Drug  Development.  Likewise,  Gene  Ontology  is  a  bio  info rmaties  tool  that 
attempts  to  unify  the  representation  of  a  gene  across  all  species  by  classifying  the  gene  in  3 
categories  whieh  inelude  1)  Cellular  Component  (CC)  2)  Molecular  Function  (MF)  or  3) 
Biological  Proecss  (BP).  Genes  that  fall  into  the  category  of  Cellular  Component  describes 
where  in  the  eell  a  gene  aets  or  what  organelle  a  gene  product  functions  in.  Moleeular  Function 
deals  with  the  elemental  activities  of  a  gene  produet  at  the  moleeular  level,  such  as  binding  or 
some  type  of  catalytie  activity.  Biological  Proeess  involves  those  genes  that  are  important  for 
the  functioning  of  eells,  tissues,  organs,  and  the  organism.  These  groupings  (i.e.,  KEGG 
pathways  and  Gene  Ontologies)  would  allow  for  a  more  complete  understanding  of  moleeular 
changes  oeeurring  during  the  stage  of  metamorphie  eompetenee  of  II.  elegans. 


B.  Genome  Sequencing  of  P.  luteoviolacea  and  C.  lytica 

PacBio  SMRT  Sequencing  of  the  P.  luteoviolacea  (HI1)  genome  resulted  in  82,296  raw  reads 
with  a  mean  read  length  of  5,344  bp,  totaling  439,749,929  nuelcotides.  Previous  genomie 
sequencing  of  this  microorganism  was  performed  by  New  Mexico  State  University  (Road 
Runner)  using  Roehe  454  sequencing  methods  (Life  Technologies)  resulting  in  133  eontigs.  To 
ereate  a  hybrid  assembly,  the  Roehe  454  sequencing  assembly  was  combined  with  the  HGAP 
assembly  using  Minimus  (Sommer  et  al.,  2007)  and  re-seaffolded  using  PaeBio’s  AHA  Program. 
The  hybrid  assembly  of  the  Roehe  454  and  PaeBio  SMRT  sequencing  resulted  in  10  scaffolds 
representing  the  5.8  Mb  genome  of  P.  luteoviolacea  (HI1). 

Similarly,  SMRT  sequencing  of  the  C.  lytica  (HI1)  genome  resulted  in  156,902  raw  reads  with  a 
mean  read  length  of  5,564bp,  totaling  873,038,51 1  nucleotides.  Generated  reads  were  then 
introduced  into  the  Hierarchical  Genome  Assembly  Proeess  (HGAP)  w'hieh  includes  assembly 
with  the  Celera  Assembler  and  assembly  polishing  with  Quiver.  The  complete  genome  of  C. 


lytica  LIM-21  was  used  as  the  reference  genome  in  the  sequence  assembly  of  C.  lytica  (HI  1 ). 
The  final  complete  genome  resulted  in  a  single  scaffold  of  3,824,196  bp,  in  comparison  to 
3,765,936  bp  of  the  published  complete  genome  of  C.  lytica  type  strain  LIM-21. 

Both  the  draft  genome  sequences  of  P.  luteoviolacea  and  the  complete  genome  of  C.  lytica  were 
submitted  to  the  RAST  (Rapid  Annotation  of  Subsystem  Technology)  server  (Aziz  et  al.,  2008) 
(Overbeek  et  al.,  2013)  which  is  a  fully  automated  on-line  server  used  to  annotate  complete  or 
nearly  complete  bacterial  and  areheal  genomes.  RAST  analysis  predieted  a  total  of  5284  coding 
sequences  for  P.  luteoviolacea  and  3340  coding  regions  for  C.  lytica.  All  coding  regions  for 
both  microorganisms  were  then  grouped  into  functional  categories  called  subsystems  (Figure  2, 
Figure  3) 


IV.  Expanded  Accomplishments 

A.  Amplification  of  Candidate  Receptor  Genes  and  In-Situ  Hybridization  of 
Competent  Larvae  of  H.  elegans. 

The  larvae  of  H.  elegans  were  grown  to  competence  (day  6)  and  harvested  on  ice,  fixed  and 
stored  in  methanol  for  in-situ  hybridization  studies.  Larvae  that  were  not  fixed  were  placed  on 
ice  and  larval  RNA  was  immediately  extracted  using  the  Qiagcn  RNcasy  Mini  Kit,  reversed 
transcribed,  and  the  cDNA  was  generated  using  the  SMARTer  RACE  cDNA  Amplification  kit 
(Clontech).  Primers  specific  to  genes  found  in  the  transcriptomc  database  and  related  to  bacterial 
recognition  sequences  were  designed  and  riboprobes  were  synthesized.  The  riboprobes  were 
then  used  in  in-situ  hybridization  experiments  to  visualize  the  anatomical  location  of  the 
transcripts  on  the  larvae  of  H.  elegans  (Figure  1). 


Figure  1.  Whole  mount  in-situ  hybridization  of  Hydroides  elegans  competent  (day  5-6)  larvae 
displaying:  (A):  Competent  larvae,  (B)  Control  (C)  Expression  of  Galactose  Binding  Lectin  (D) 
Expression  of  Peptidoglycan  Recognition  Protein  (E)  Expression  ofNF-KappaB 

B.  Differential  Expression  of  Genes  Unique  to  the  Competent  Larvae  of  H.  elegans 

The  larvae  of  H.  elegans  was  grown  to  precompetence  (Day  3),  competence  (day  6)  and  Adult 
stages  and  harvested  on  iee.  The  larval  and  adult  RNA  was  then  extracted  using  the  Qiagcn 
RNeasy  Mini  Kit.  The  quality  of  RNA  was  then  determined  using  the  Agilent  Bioanalyzer  and 


submitted  to  the  Huntsman  Cancer  Institute  at  the  University  of  Utah  for  transcriptome 
sequencing.  After  the  alignment  and  assembly  of  the  reads,  all  sequences  were  submitted  to 
Blast  and  compared  against  databases  at  NCBI,  Swiss-Prot,  TRembl,  and  NR  All  the  identified 
protein  sequences  that  were  differentially  expressed  in  only  competent  larvae  were  then 
determined  using  EBSeq  and  called  Pattern  3  (Table  1).  Pattern  3  gene  sequences  were  then 
analyzed  using  DAVID  and  classified  into  KEGG  Pathways  (Table  2)  and  Gene  Ontologies 
(GO)  (Table  3).  Table  2  shows  the  most  significant  biological  pathways  (p<0.05)  involved  at  the 
stage  of  metamorphic  competence.  Likewise,  Table  3  shows  the  most  significant  processes 
involved  in  the  GO  categories  of  Cellular  Component,  Molecular  Function,  and  Biological 
Process.  Overall,  KEGG  Pathway  analysis  shows  that  at  the  stage  of  metamorphic  competence 
there  is  an  upregulation  of  genes  that  are  involved  in  processing  of  environmental  information  by 
way  of  genes  involved  in  Neuroactive  Ligand  Receptor  and  Notch  Signaling  Pathways. 

Likewise,  on  an  organismal  level,  there  is  an  up-regulation  of  genes  involved  in  the  development 
(Dorso-Ventral  Axis  Formation),  circulatory  system  (Vascular  Smooth  Muscle),  and  the  immune 
system  (Complement  and  Coagulation  Pathway).  Gene  Ontology  results  also  show  that  genes 
that  are  significantly  expressed  during  metamorphic  competence  involve  receptor  signaling 
activity  (Biological  Process  -  GPCR,  Molecular  Function  -  Neurotransmitter  Activity). 

Table  1.  Various  pattern  designations  for  the  differential  expression  of  genes  of 
precompetent,  competent  and  adult  stages  of  H.  elegans. 


Pattern 

Gene  Expression 

Precompetent 

Competent 

Adult 

1 

Same  for  all  stages 

1 

1 

1 

2 

Unique  to  precompetent  stage,  common  to 
competent  and  adult  stages 

2 

1 

1 

3 

Unique  to  competent  stage,  common  to 
precompetent  and  adult  stages 

1 

2 

1 

4 

Unique  to  adult  stage,  common  to  precompetent 
and  competent  stages 

2 

2 

1 

5 

Different  in  all  stages 

3 

2 

1 

Table  2.  KEGG  Pathways  involved  at  the  stage  of  competence  for  H.  elegans 


KEGG  Pathway 

Category 

P-value 

Neuroactive  Ligand  Receptor  Interaction 

Environmental 
Information  Processing 

9  x  10'5 

Vascular  Smooth  Muscle 

Organismal 

6  x  10'J 

Notch  Signaling  Pathway 

Environmental 
Information  Processing 

2.4  x  10'" 

Complement  and  Coagulation  Cascade 

Organismal 

2.4  x  10'" 

Proteosome 

Genetic  Information 
Processing 

2.5  x  10‘" 

Dorso-Ventral  Axis  Formation 

Organismal 

3.8  x  10'" 

Chondroitin  Sulfate  Biosynthesis 

Metabolism 

4.6  x  10’" 

Table  3.  Overview  of  Gene  Ontology  Classification  for  Competent  Larvae 


Gene  Ontology 

Process 

P-value 

Cellular  Process 

Integral  to  Plasma  Membrane 

2.4  x  10’15 

Extracellular  Region 

8.8  x  KT* 

Synapse  Region 

5.4  x  fir7 

Molecular  Function 

Calcium  Ion  Binding  Activity 

4.4  x  10'17 

Guanylate  Cyclase  Activity 

1.5  x  10’12 

Neurotransmitter  Receptor  Activity 

6.9  x  10'iU 

Biological  Process 

G-Protein  Coupled  Receptor  Signaling  Pathway 

6.8  x  10'* 

cGMP  Biosynthesis 

5.7  x  10'u 

Carotenoid  metabolism 

1.7  x  10'2 

C.  RAST  Analysis  of  P.  luteoviolacea  and  C.  lytica  Genomes 

A  recent  study  on  P.  luteoviolacea  has  found  that  the  inductive  gene  products  responsible  for  the 
settlement  and  metamorphosis  of  H.  elegans  are  components  of  a  phage  tail-like  bacteriocin 
(Shikuma  et  al.,  2014).  RAST  analysis  of  the  P.  luteoviolacea  genome  (Figure  2)  shows  that 
there  are  14  genomic  elements  within  the  genome  of  P.  luteoviolacea  that  are  found  in  the 
category  Phages,  Prophages,  Transposable  Elements  and  Plasmids.  Most  of  these  sequences 
encode  for  phage  structural  components  such  as  the  sheath,  tube  and  base  plate  of  a  phage  tail. 
This  is  in  contrast  to  the  C.  lytica,  the  second  inductive  strain,  which  has  2  elements  related  to 
Phages,  Prophages,  Transposable  element  and  Plasmids  contained  within  its  genome  as 
determined  by  RAST  (Figure  3).  The  two  genomic  elements  of  C.  lytica  encode  for  a  phage-tail 
fiber  protein  and  a  Listeria  Pathogenicity  Island,  but  no  other  phage-tail  structural  component  is 
identified  (i.e.,  phage  sheath,  tube  or  base  plate).  Thus,  it  is  appears  that  the  cue  responsible  the 
settlement  and  metamorphosis  of  H.  elegans  produced  by  C.  lytica  involves  a  different 
mechanism  and  requires  further  analysis  and  research  in  the  future. 


Subsystem  Category  Distribution 


Subsystem  Feature  Counts 

g  ■  Cofactors,  Vitamins,  Prosthetic  Groups,  Pigments  (224) 
gfl  CeK  Wall  and  Capsule  (143) 
g  ■  Virulence,  Disease  and  Defense  (108) 

Potassium  metaboism  (19) 

Photosynthesis  (0) 

©  Miscelaneous  (21) 

0  ■  Phages,  Prophages,  Transposable  elements,  Plasmids  (14) 
£)■  Membrane  Transport  (148) 

SB  Iron  acquisition  and  metabolism  (37) 

(3  RNA  Metabolism  (186) 
g  ■  Nucleosides  and  Nucleotides  (88) 

Q  ■  Protein  Metabolism  (270) 

O  a  Cei  Division  and  Ceil  Cycle  (29) 
q  Motiity  and  Chemotaxis  (134) 
a  ■  Regulation  and  Cell  signaling  (79) 

0  Secondary  Metabolism  (10) 

QH  DN A  Metabolism  (115) 

@ ■  Regulons  (0) 

iB  ■  Fatty  Acids,  Lipids,  and  Isoprenoids  (112) 

0B  Nitrogen  Metabolism  (13) 

0fl  Dormancy  and  Sporulation  (4) 
gH  Respiration  (91) 

0  Stress  Response  (149) 
a  Metabolism  of  Aromatic  Compounds  (18) 
a  Amino  Acids  and  Derivatives  (379) 
g  Sulfur  Metabolism  (30) 

£!■  Phosphorus  Metabolism  (41) 

3 ■  Carbohydrates  (259) 


Figure  2:  RAST  analysis  of  Pseudoaheromottas  luteoviolacea 


Subsystem  Category  Distribution 


Subsystem  Feature  Counts 

Cofactors,  Vitamins,  Prosthetic  Groups,  Pigments  (182) 
Cel  Wall  and  Capsule  (92) 

Virulence,  Disease  and  Defense  (53) 

Potassium  metaoolism  (17) 

Photosynthesis  (0) 

Miscellaneous  (21) 

Phages,  Prophages,  Transposable  elements.  Plasmids  (2) 
Membrane  Transport  (79) 

Iron  acquisition  and  metabolism  (4) 

RNA  Metabolism  (125) 

Nucleosides  and  Nucleotides  (56) 

Protem  Metabolism  (148) 

Cell  Division  and  Cell  Cycle  (17) 

Motility  and  Chemotaxvs  (O) 

Regulation  and  Cell  signaling  (17) 

Secondary  Metabolism  (7) 

DNA  Metabolism  (70) 

Regulons  (5) 

Fatty  Acids,  Upids,  and  Isoprenoids  (55) 

Nitrogen  Metabolism  (17) 

Dormancy  and  Sporulation  (3) 

Respiration  (48) 

Stress  Response  (59) 

Metabolism  of  Aromatic  Compounds  (19) 

Amino  Acids  and  Derivatives  (276) 

Sulfur  Metabolism  (33) 

Phosphorus  Metabolism  (18) 

Carbohydrates  (213) 


Figure  3.  RAST  Analysis  of  Cellulophaga  lytica 
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